arXiv: 1501.07327vl [physics.atom-ph] 29 Jan 2015 


Spectra of barium, radium, and element 120; application of the combined correlation 
potential, singles-donbles, and confignration interaction ab initio method 

J. S. M. Ginges and V. A. Dzuba 

School of Physics, University of New South Wales, Sydney NSW 2052, Australia 

(Dated: January 30, 2015) 

We apply a version of the recently developed approach combining the correlation potential, lin¬ 
earized singles-doubles coupled-cluster, and the configuration interaction methods to the spectra of 
the heavy alkaline earths barium, radium, and element 120. Quantum electrodynamics radiative 
corrections are included. We have found unprecedented agreement between ab initio theory and 
experiment for the spectra of barium and radium, and we make accurate predictions for missing and 
unreliable data for all three atoms. 


I. INTRODUCTION 

Preparatory work with the radium atom towards mea¬ 
surements of fundamental symmetries violations is in 
progress at Argonne National Laboratory (USA) [ij and 
Kernfysisch Versneller Instituut (The Netherlands) Q. 
Studies of the parity and time-reversal violating atomic 
electric dipole moment (EDM) and atomic parity viola¬ 
tion (APV) are particularly attractive in radium due to 
orders of magnitude enhancement of the effects, arising 
from both nuclear and electronic mechanisms: the pres¬ 
ence of octupole deformation of the nucleus may lead to 
several-hundred times nuclear enhancement of the EDM 
in radium in the electronic ground state compared to 
mercury (see, e.g., Si) , for which the best limit on an 
atomic EDM has been placed i ; the presence of anoma¬ 
lously close electronic levels of opposite parity may lead 
to orders of magnitude electronic enhancement of EDM 
effects in metastable states and APV effects in certain 
transitions Si- 

High-precision atomic structure calculations will be re¬ 
quired for interoretation of the APV and atomic EDM 
measurements (8|. Moreover, the measured low-lying ex¬ 
citation energies for radium are incomplete, and it is im¬ 
portant to have high-precision predictions for these levels. 
There are already a number of theoretical works devoted 
to studies of fundamental symmetries violations in ra¬ 
dium BIMl and to the radium spectra, lifetimes, and 
hyperfine structure [Tsl - I^ . There are also recent mea¬ 
surements of a few transition frequencies and lifetimes 
(23 - 1^ . In this work we use a recently-developed ap¬ 
proach that combines the correlation potential, singles- 
doubles coup led-cluster, and configuration interaction 
methods to the radium spectra. We have found un¬ 
precedented agreement between theory and experiment 
for most levels using this ab initio method. 

The heavier electronic homologue of radium is element 
120. Efforts to synthesize this element are underway at 
GSI (Germany), JINR (Russia), and RIKEN (Japan) 
(see, e.g., [1^). This element lies within the predicted 
island of stability, a region of increased stability against 
nuclear decay close to the next doubly magic shell clo¬ 
sures Z=114, 120, or 126 and N=172 or 184, depend¬ 
ing on the model [l^. If Z=120 is the next closed pro¬ 


ton shell, significantly increased stability is expected for 
this element. Enhanced stability could make atomic and 
chemical studies of this element possible. Already, chem¬ 
ical studies involving the superheavy element Sb (Z=106) 
have proved successful [l^. The current calculations for 
E120 extend the spectral range considered in previous 
works [29l - l^ . 

We also perform calculations for the lighter electronic 
homologue barium, as this can be used to gauge the ac¬ 
curacy for the heavier elements. The method used in this 
work has alread y b een applied to the low-lying levels of 
barium in Ref. |25j |. Here we extend the application to 
the excitation energies of the lowest 31 levels. 

II. METHOD OF CALCULATION 

We use an approach that is based on the combination 
of several different many-body methods: the correlation 
potential (CP), linearized singles-doubles coupled-cluster 
(SD), and configuration interaction (Cl) methods. It may 
be referred to as the CP-I-SD-I-CI method. This method 
was recently developed by Dzuba in Ref. [1^ and is 
similar to the combined SD-I-CI method developed by 
Safronova et al. [s^ . 

The CP-l-SD-l-CI and SD-I-CI methods are essentially 
based on the method combining many-body perturbation 
theory and the configuration interaction (MBPT-I-CI) 
[s^ . which has proven to be one of the most compu¬ 
tationally efficient and accurate approaches for calcula¬ 
tions involving heavy atoms with several valence elec¬ 
trons. These methods differ in their treatment of the 
valence-core correlations and the screened Coulomb in¬ 
teraction (valence-core-valence correlations), as we ex¬ 
plain below. 

Barium, radium, and element 120 have two valence 
electrons above a closed electronic core. We wish to find 
the eigenvalues and eigenvectors of the effective Hamil¬ 
tonian for the two valence electrons, 

= hi{ri)+hi{r2) + h2{ri,r2) , ( 1 ) 

by diagonalizing with respect to wavefunctions con¬ 
structed from linear combinations of two-electron Slater 
determinants; the Slater determinants are formed from 
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TABLE I. Factors for Ba used to mimic higher-order screen¬ 
ing for exchange diagrams, fk = ,.|?;) and 

fk’^ = (^|Sdi’r,fe|wy(’^|Sdn,fel^); ^ i® ^he multipolarity of the 
Coulomb interaction. 


k 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

fk 

0.73 

0.61 

0.83 

0.91 

0.99 

1.05 

1.10 

1.00 

1.00 

1.00 

rhp 

Jk 

0.86 

0.73 

0.97 

1.02 

1.07 

1.11 

1.15 

1.00 

1.00 

1.00 


relativistic Hartree-Fock orbitals found in the core of 
N — 2 electrons. Here hi contains all one-electron terms 
of the many-electron problem and h 2 contains all two- 
electron terms. The one-electron terms are 

hi= ca-p + {P - l)mc^ -h t4uc + ’ (2) 

where a. and P are Dirac matrices, Hnuc is the nuclear 
potential (we use the two-parameter Fermi distribution 
for the nuclear density), is the self-consistent 

Hartree-Fock potential of the TV — 2 electrons making 
up the electronic core, and Ei is the correlation poten¬ 
tial which accounts for the correlations between a single 
valence electron and the core. For E120, the choice of 
the nuclear parameters is important; for our Fermi dis¬ 
tribution, we use a nuclear radius of r = 8.0 fm and nu¬ 
clear thickness t = 2.0 fm, corresponding to a root-mean- 
square radius rrms ~ 6.42 fm. These are the parameters 
we used in our previous calculation [2^, and we make 
this choice again for easy comparison with that work. 

The two-electron terms ^2 consist of the Coulomb in¬ 
teraction and the two-electron correlation operator E 2 
(the screened Coulomb interaction), 

^2(r-i,r2) =-r + S2(ri,r2). (3) 

Id - d| 

In the MBPT-I-CI method, both Ei and E 2 are cal¬ 
culated using MBPT; they are often calculated in only 
second-order in the Coulomb interaction. We have 
demonstrated before, that using the all-order correlation 
potential E“ in place of the second-order correlation po¬ 
tential leads to significantly improved accuracy [l^ . 
We refer to this latter approach, with the all-order 
as the CP-I-CI method. In the SD-I-CI method [s^, both 
El and E 2 are found from all-order linearized singles- 
doubles coupled-cluster equations. 

In the current CP-l-SD-fCI approach, we calculate 
the all-order E“ using the Feynman diagram technique, 
while the screening of the Coulomb interaction E 2 is 
found from the SD method. In this way, we include im¬ 
portant classes of many-body diagrams to all-orders in 
the Coulomb interaction, both for Ei and E 2 . It was 
shown in Ref. that high accuracy in excitation ener¬ 
gies can be obtained with this choice. 

In the calculation of Ei, the Feynman diagram tech¬ 
nique is used to include two classes of diagrams to all 
orders in the Coulomb interaction for the direct part: 


the electron-electron Coulomb screening and the hole- 
particle interaction inside the internal loops [s^. The ex¬ 
change part is evaluated in second-order in the Coulomb 
interaction by calculation of the corresponding Goldstone 
diagrams. Screening is taken into account in a simplified 
way, by multiplying the Coulomb integrals by factors 
found from the direct diagrams (more on this below). 
There is another series of diagrams, referred to as “lad¬ 
der diagrams”, that are calculated using singles-doubles- 
type equations [s^ ; the corresponding corrections to the 
valence energies are added to the Cl matrix. 

Screening factors for the exchange diagrams are found 
by taking the ratio of the expectation value for the di¬ 
rect part of (with the dominant all-orders electron- 

electron screening considered) to the expectation value 
for the direct part of E) ' for each multipolarity fc, that 
is, fk = The corresponding fac¬ 

tors for Ba are shown in the first row of Table HI Later, 
as a test of our accuracy, we consider the effect of using 
a different set of screening factors found by consid¬ 
ering both the dominant electron-electron screening and 
the hole-particle interaction series of diagrams; these fac¬ 
tors are listed in the second row of Table HI 

Note that we used slightly different factors for fk in 
different approximations (with or without Breit or QED 
corrections) as well as slightly different factors for differ¬ 
ent atoms (Ba, Ra, E120). They differ by a few percent 
at most, and as their precise values are not of general 
interest, we do not present them here. 

We refer the reader to Ref. [ 2 ^ for details regarding 
the calculation of E 2 . 

In this work we quantify the corrections associated 
with inclusion of the Breit interaction and the quantum 
electrodynamics radiative corrections. The Breit interac¬ 
tion is considered at the relativistic Hartree-Fock (RHF) 
level. The radiative corrections are taken into account 
through the addition of a local radiative potential [s^ to 
the nuclear potential. This radiative potential method 
provided a breakthough in the ability for radiative cor¬ 
rections to be included into the many-body problem in 
heavy atoms. We refer the reader to Ref. for details 
about this method. The radiative potential method has 
since been implemented in a number of works, including 
(with a minor modification to one of the fitting factors 
in the radiative potential) Refs. [13, where the ra¬ 
diative potential was added to the atomic many-body 
package GRASP [H, . 

The upper and lower radial components of the rela¬ 
tivistic Hartree-Fock orbitals in the V^~‘^ potential are 
expanded in a basis of B-splines of order fc = 9 [HI . We 
use 40 B-splines for Ba and Ra and 50 B-splines for E120 
confined to a cavity of radius 40 a.u. We perform the 
calculations for the exchange part of the correlation po¬ 
tential, the ladder diagrams, and E 2 using the lowest 30 
states in each partial wave up to / = 6 for the intermedi¬ 
ate excitations. We include the correlation potential Ei 
for the valence levels up to T = 4 and E 2 for the valence 
levels up to T = 3. At the Cl stage of the calculations. 
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TABLE II. Ionization energies for the lowest states in each 
wave to I — 3 for Ba^; units cm“^. 


State 

RHF 

-f 

Lad. 

Breit 

QED 

Total 

Exp.^ 

A 

6Si/2 

75340 

80808 

-156 

-22 

-45 

80585 

80686 

-101 

5^3/2 

68139 

76466 

-763 

34 

36 

75773 

75812 

-39 

5^5/2 

67665 

75622 

-765 

57 

32 

74946 

75011 

-65 

6pi/2 

57266 

60591 

-128 

-36 

5 

60432 

60425 

7 

6P3/2 

55873 

58865 

-118 

-16 

3 

58734 

58734 

0 

4 / 5/2 

28213 

32180 

-32 

117 

32 

32297 

32428 

-131 

4 / 7/2 

28222 

31989 

-32 

126 

28 

32112 

32203 

-91 


NIST data, Ref. [H. 


the basis for Ba and Ra (and E120) consists of the lowest 
14 (18) valence states in each partial wave up to 1 = 4. 

To quantify the size of the Breit and QED corrections, 
we have carried out three runs of the calculations, for 
both the ions and the neutral divalent atoms: with nei¬ 
ther Breit nor QED, with Breit and without QED, and 
with both Breit and QED. As the basis is determined 
by the Hartree-Eock orbitals, and Breit and QED correc¬ 
tions enter at this level, we have used a different basis set 
for each of these runs. For each run, one basis is used for 
all aspects of the many-body problem: CP, SD, and CL 


III. RESULTS 

Our calculations begin for the ions Ba’*', Ra’*', and 
E120'''. The correlation potential Ei that describes the 
valence-core correlations is the same for the monovalent 
ions above as for the divalent atoms Ba, Ra, E120 when 
we perform the RHF calculations in the field of the N — 2 
electrons of the core. Therefore, the quality of the spec¬ 
tra for the ions is a good indication of the quality of the 
correlation potential Ei (T^ . 

A. Ions 

Results for the ionization energies of the lowest partial 
waves up to Z = 3 for the ions Ba’*', Ra’*', and EI20'*' 
are presented in Tables miiiii and IlYl respectively. We 
separate the contributions arising from inclusion of the 
correlation potential (without ladder diagrams), lad¬ 
der diagrams “Lad.”, the Breit interaction, and QED ra¬ 
diative corrections. 

Absolute differences of our final results from experi¬ 
ment are listed in the last column under “A”. These dif¬ 
ferences for Ba"*" and Ra’*' are very small, on the order of 
10—100cm“^, with the largest difference |A| ss 160 cm“^. 
We expect a similar level of agreement for E120“''. 

We can see from a comparison of the final results “To¬ 
tal” in Tables |TT1 HID IIVI the effect of the relativistic con¬ 
traction of the s and Pi /2 orbitals in the heavier homo- 
logues, which in turn screen and push out the orbitals 
with higher angular momenta. 


Inclusion of the ladder diagrams is very important for 
reaching good agreement for the d levels. The d-orbitals 
are very low-lying in the ions and play a big role in the 
low-lying spectra of the neutral divalent atoms. The size 
of the ladder diagram contribution decreases as we go 
from Ba"*" to E120''', as the d orbitals are pushed out, 
while the contributions for s and p levels increase. 

The contributions from Breit and QED are roughly of 
the same magnitude, and generally increase with higher 
Z. The QED corrections exceed 100 cm“^ for the s lev¬ 
els for E120+, while they are negligible for the p waves. 
The largest Breit corrections we see are for the / levels, 
almost reaching 200 cm“^. Interestingly, these /-wave 
corrections are mostly determined by many-body effects 
through the inclusion of and they are sensitive to 
the choice of the correlation potential. 


TABLE III. Ionization energies for the lowest states in each 
wave to Z = 3 for Ra"^; units cm“^. 


State 

RHF 

-f 

Lad. 

Breit 

QED 

Total 

Exp.^ 

A 

7Si/2 

75898 

82010 

-219 

-21 

-87 

81684 

81842 

-158 

6^3/2 

62356 

70242 

-620 

46 

38 

69707 

69758 

-51 

6^5/2 

61592 

68518 

-643 

67 

31 

67973 

68099 

-126 

7Pl/2 

56878 

60739 

-182 

-59 

2 

60499 

60491 

8 

7P3/2 

52906 

55771 

-140 

-18 

-2 

55611 

55633 

-22 

5 / 5/2 

28660 

32768 

-74 

109 

32 

32835 

32854 

-19 

5 / 7/2 

28705 

32542 

-76 

107 

27 

32600 

32570 

30 


^ NIST data. Ref. O. 


It is worthwhile pointing out that the radiative poten¬ 
tial [1^ was found by fitting to the s and p levels for 
the hydrogen-like ions; it was not fitted to higher waves. 
The radiative QED corrections for the d and / waves, 
however, largely arise due to account of many-body ef¬ 
fects, namely core relaxation and the correlation poten¬ 
tial. Core relaxation corresponds to alteration of the elec¬ 
tronic core due to self-consistent solution of the RHF 
equations with the radiative potential included. The 
largest part of the core relaxation correction to the va¬ 
lence level shift occurs as a result of radiative corrections 
to s-orbitals in the core, which are well-fitted. 

We have previously calculated the spectra of these ions 
for s and p levels, including estimates of Breit and radia¬ 
tive corrections ; more recently the spectra was cal¬ 
culated in Ref. for the lowest s, p, and d levels, this 
time with the contribution from ladder diagrams, also. 

B. Divalent atoms 

In Tables IVlIVIl and I VIII we present our results for the 
ionization potentials (removal of one s electron, IPl, and 
removal of both s electrons, IPI-I-IP2) and excitation en¬ 
ergies for the lowest 31 levels for Ba, lowest 40 levels for 
Ra, and lowest 25 levels for E120. Results presented un¬ 
der the column “CP"*-|-SD-|-Cr’ mean that the all-orders 
correlation potential is included (though ladder di- 
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TABLE IV. Ionization energies for the lowest states in each 
wave to I = 3 for E120’*'; units cm“^. 


State 

RHF 

± 

Lad. 

Breit 

QED 

Total 

00 

83168 

90105 

-520 

-108 

-129 

89349 

8 pi /2 

60027 

65475 

-379 

-141 

-9 

64946 

7d3/2 

56620 

64841 

-593 

51 

58 

64357 

7d5/2 

56413 

62684 

-625 

46 

46 

62151 

8P3/2 

49295 

52017 

-179 

-18 

-7 

51812 

6 / 5/2 

29734 

36907 

-211 

186 

97 

36978 

6 / 7/2 

29909 

36252 

-219 

169 

77 

36279 


agrams are not taken into account; the superscript “nl” 
is short for “no ladder”) and that the all-orders E 2 is 
included (calculated using the SD method). Contribu¬ 
tions from ladder diagrams, Breit, and QED radiative 
corrections appear in the following columns. 

Our final results are presented under the column “To¬ 
tal” , and the column “Aexp” gives the deviation of these 
results from experiment, Aexp = Total — Exp. Results 
of other calculations are given in the tables, also. In 
particular, we have presented the most precise ab initio 
calculations available. These include different versions of 
coupled-cluster [H, HO, HH Sll , CP-I-CI [ 2 ^, and SD-I-CI 
[S^. For the higher levels, where there is only limited 
data available, we also present for Ba and Ra the results 
of semi-empirical CP-I-CI calculations [^ . 


1. Barium 

For Ba, it is seen that the ladder diagrams give a con¬ 
tribution to the excitation energies of about 500 cm“^ for 
configurations containing a single 5d orbital, that is for 
6s5d and 5d6p (excluding the very highest level). For the 
bcP configuration, the contribution amounts to roughly 
double this, ~ 1000 cm“^. The ladder diagrams give a 
small correction to the other levels. 

The Breit and radiative corrections for Ba are roughly 
of the same magnitude, ranging between about 10 — 
120 cm“^ for Breit and 15 — 110 cm“^ for the radiative 
corrections for the levels considered. For the higher levels 
the radiative corrections dominate. 

We see that the deviation of our final results “Total” 
from experiment, Aexp, ranges between 10 — 200 cm“^, 
with the exception of the larger deviation for bcP ^Po,i of 
about 400 cm“^ and 500 cm“^ and the very large devia¬ 
tion for of about 1100 cm“^. 

The singlet state has not been measured and 

is absent in the NIST data [l^. Its position was pre- 
dicted in the early experimental work of Palenius to 
be 24300 ±300 cm“^. It has subsequently been calculated 
in the works [l^, I47l - l4^ ; there are other calculations of 
the barium spectra where this level has been missed. We 
present in the table, alongside our own result, the initial 
estimat e 14611 and the value from the most precise calcu¬ 
lations [ 3 . Our calculations give the value 25205 cm“^ 


for this level; looking at the results for other terms with 
the configuration 5d^, we expect that the result for this 
level could be underestimated, possibly by as much as 
100- 300cm-k 

Our result for 5d^ disagrees with the experimental 
value by k, 1100cm“^. This very large difference is well 
outside the deviations we see for the other levels. This 
level does not appear in the extensive spectra calculations 
of Landau et al. [l^ . We know of only one other calcula¬ 
tion of this level, carried out in the CP±CI method with 
empirical fitting (2^, with the result 26034 cm“^; this is 
about 700 cm“^ less than experiment, and well above the 
estimated error in that work. 

We should note, however, that the largest deviations 
we see are for terms belonging to 5d^, and the unmea¬ 
sured level bd? and the level bcP ^Sq belong to this 
configuration. It is possible that our calculations do not 
do an adequate job for these low-lying (P states. We note 
further the large difference between our result and that 
of Ref. [3 for the energy of the level b(P ^64 . 

The results of the current work are in significantly bet¬ 
ter agreement with experiment over the considered spec¬ 
tral range compared to the ab initio calculations per¬ 
formed using CP±CI [ 2 ^, SD±CI 1^. and a version of 
Fock-space coupled-cluster (FSCC) [^. There is only 
one other work where comparable accuracy was reached, 
using the intermediate Hamiltonian Fock-space coupled- 
cluster (IHFSCC) method [l^; however, in that method 
the accuracy deteriorates for the higher levels. 

We should note that neither Breit nor QED corrections 
were included in the results of Ref. [ 2 ^ , while QED cor¬ 
rections were not included in those of Refs. [H, 1^.14511. 
The QED correction to IPl was calculated in the work 
[ 3 ^ using the radiative potential [s^, and the value 
— 19cm“^ was obtained - the same result we have ob¬ 
tained in this work. 


2. Radium 

The results for radium are presented in Table |Vll It is 
seen from the table that there are a number of gaps in the 
experimental data, and accurate theoretical predictions 
of the missing data are important. We don’t agree with 
all configuration designations used in the experimental 
data compiled by NIST [i^, and we predict that the 
energy of one of the states {7s8p ^Pi) is significantly 
lower than that given in the data tables, as explained 
later. 

Compared to Ba, the ladder contributions for Ra are 
smaller for the terms involving d orbitals in the dominant 
configurations, while the contributions to terms from the 
ss and sp configurations are larger. The Breit contribu¬ 
tions to the excitation energies range from 10 — 81 cm“^ 
for the levels considered, while we saw a significantly 
larger contribution arising from the QED radiative cor¬ 
rections, 37 — 180 cm“^. 

Unlike the case with Ba, for Ra the term designations 
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TABLE V. Ionization potentials (removal of one electron, IPl, and both electrons, IP1+IP2) and excitation energies for 
Ba. CP"*+SD+CI excludes ladder diagrams, Breit, and QED; these corrections are given in subsequent columns. Aexp = 
Total — Exp. Results of other calculations are presented in the hnal column. Units: cm“^. 


Conf. 

Term 

J 

Cpnl 

+SD+CI 

Ladder 

Breit 

QED 

Total 

Exp.^ 

A Exp 

Other 

Tpi 


0 

42097 

-173 

122 

-19 

42027 

42035 

-8 

42444^, 42120^ 

IP1+IP2 6s^ 


0 

122905 

-329 

100 

-64 

122612 

122721 

-109 

122954^, 123363^, 123869^ 

6s5d 


1 

8460 

575 

-70 

-60 

8905 

9034 

-129 

9117^, 8687^, 9249^, 9077^ 



2 

8704 

576 

-68 

-59 

9153 

9216 

-63 

9296^, 8875^, 944ffi, 9369^ 



3 

9066 

576 

-61 

-57 

9524 

9597 

-73 

9677^, 9279^, 9840^, 9830^ 

6s5d 


2 

10996 

589 

-62 

-66 

11457 

11395 

62 

11426^, 1108+, 11727^, 1187T 

6 s6p 

3^0 

0 

12351 

-33 

27 

-35 

12310 

12266 

44 

12357^,12099^, 12556^, 12668^ 



1 

12706 

-33 

26 

-35 

12664 

12637 

27 

12728^, 12474_^, 12919^, 1294^ 



2 

13603 

-41 

25 

-32 

13555 

13515 

40 

13610^, 13365^, 13819^, 13449J 

6 s6p 

IpO 

1 

17887 

107 

11 

-41 

17964 

18060 

-96 

18170^, 17943^, 18292^, 20077J 

5(f 

3p 

2 

19920 

1142 

-119 

-110 

20833 

20934 

-101 

21017^ 



3 

20142 

1146 

-115 

-108 

21065 

21250 

-185 

21338^ 



4 

20636 

1145 

-109 

-106 

21566 

21624 

-58 

21714^ 

5d6p 

3p° 

2 

21545 

545 

-46 

-96 

21948 

22065 

-117 

22238^ 



3 

22436 

538 

-42 

-93 

22839 

22947 

-108 

23116^ 



4 

23286 

534 

-39 

-91 

23690 

23757 

-67 

23950^ 

bcP 


2 

22182 

996 

-86 

-103 

22989 

23062 

-73 

23077^ 

5d6p 


2 

22607 

543 

-48 

-94 

23008 

23074 

-66 

23289^ 


3p 

0 

21985 

1001 

-88 

-104 

22794 

23209 

-415 

23213^ 



1 

22118 

1022 

-90 

-104 

22946 

23480 

-534 

23500^ 



2 

22897 

1003 

-81 

-100 

23719 

23919 

-200 

23950^ 

5d6p 


1 

23635 

541 

-61 

-98 

24017 

24192 

-175 

24474^ 



2 

23992 

545 

-58 

-97 

24382 

24532 

-150 

24817^ 



3 

24460 

543 

-56 

-95 

24852 

24980 

-128 

25264,^ 

bd^ 


4 

24328 

1069 

-89 

-103 

25205 



24300(300)1, 246841. 

5d6p 

3pO 

0 

25115 

507 

-54 

-94 

25474 

25642 

-168 

258861. 



1 

25183 

503 

-51 

-92 

25543 

25704 

-161 

259471. 



2 

25466 

512 

-46 

-90 

25842 

25957 

-115 

262031. 

6s7s 


1 

26253 

-135 

28 

-15 

26131 

26160 

-29 

265731. 

bd? 


0 

25144 

612 

-28 

-70 

25658 

26757 

-1099 

260341. 

bdQp 


3 

26343 

498 

-56 

-94 

26691 

26816 

-125 

269681. 

&s7s 


0 

28192 

55 

10 

-37 

28220 

28230 

-10 

285831. 

bd6p 

1 

1 

28351 

210 

-15 

-59 

28487 

28554 

-67 

287881. 


NIST data, Re fJ^ . 

^ IHFSCC, Ref. EpP 

CP+CI, Ref. [2^. IPl found in combination with the calculation for the ion, Ref. [4,'lll . 
SD+CI, Ref H. 

® FSCC, Ref. [HP 

^ Estimated in the experimental work of Palenius, Ref. 1461 . 

® Semi-empirical CP+CI, Ref. ISoll . 


change in some cases when we go from the approximation 
without ladder diagrams, CP"^+SD+CI, to that with 
ladder diagrams. This means that it is more difficult to 
track the changes in energies associated with the ladder 
corrections, as the terms themselves may differ in these 
approximations. 

The deviations from experiment are generally smaller 
for Ra than for Ba, with the deviations well under 
100 cm“^ for most levels; only a few levels deviate more 
than 100 cm“^, with the maximum deviation 191 cm“^ 
for the singlet state 7s6d^D2 (with the exception of 
7s8p ^Pi, which we address below). 


In the final column of Table IVll we list the results of 
other calculations. For the lower levels, we present the 
results of ab initio calculations, namely from the CP+CI 
method [2^ and the extended intermediate Hamiltonian 
Fock-space coupled-cluster (XIHFSCC) method [Slj. For 
the higher levels, where data from these methods is un¬ 
available, we present the results of the IHFSCC method 
[Tsil and a semi-empirical CP+CI calculation (2^. Both 
of these calculations miss some of the higher levels that 
we see. Note that QED corrections to the excitation en¬ 
emies were not taken into account in these other works 
[isl [ 2 ^, [Slj , while Breit corrections were also omitted in 
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TABLE VI. Ionization potentials (removal of one electron, IPl, and both electrons, IP1+IP2) and excitation energies for 
Ra. CP“*+SD+CI excludes ladder diagrams, Breit, and QED; these corrections are given in subsequent columns. Asxp = 
Total — Exp. Results of other calculations are presented in the final column. Units: cm“^. Asterisks (*) identify where 
configurations have been modified from those in the NIST data [d^] or where calculated energies in the final column have been 
reassigned to different terms. 


Conf. 

Term 

J 

Cpnl 

-bSD-bCI 

Ladder 

Breit 

QED 

Total 

Exp. 2 

Aexp 

Other 

1 p1 7s-‘ 


0 

42680 

-227 

271 

-45 

42679 

42573 

102 

42622^, 42562^ 

IPl-bIP2 

Is 

0 

124690 

-446 

250 

-132 

124363 

124416 

-53 

124656^, 124642^ 

7s7p 

3pO 

0 

13173 

-43 

70 

-64 

13136 

13078 

58 

12916^, 13093^ 



1 

14080 

-47 

71 

-64 

14040 

13999 

41 

13844_^, 14017^ 



2 

16828 

-84 

78 

-60 

16762 

16689 

73 

16566^, 16675^ 

7s6d 


1 

13411 

380 

-46 

-92 

13653 

13716 

-63 

13622_^, 1402R 



2 

13771 

384 

-39 

-91 

14025 

13994 

31 

13902_^, 14292,^ 



3 

14440 

401 

-20 

-85 

14736 

14707 

29 

14645^, 14989^ 

7s6d 


2 

16996 

385 

-12 

-97 

17272 

17081 

191 

17004^, 17376^ 

7s7p 


1 

20632 

-34 

69 

-61 

20606 

20716 

-110 

20667^, 20792,^ 

7s8s 


1 

26887 

-181 

73 

-37 

26742 

26754 

-12 

26665^, 26762^ 

7s8s 

is 

0 

27910 

-119 

70 

-45 

27816 



27768^, 28248^ 

6d7p 

3pO 

2 

27803 

337 

26 

-161 

28005 

28038 

-33 

27991^, 28328^ 



3 

29924 

328 

46 

-155 

30143 

30118 

25 

30067^, 30388^ 



4 

32248 

303 

59 

-151 

32459 

32368 

91 

32363^, 32603^ 

&d? 

3p 

2 

28964 

738 

-67 

-180 

29455 



29731^, 29610^ 



3 

29648 

776 

-56 

-179 

30189 



30464_^, 30404,^ 



4 

30455 

787 

-39 

-175 

31028 



31172^, 31114^ 

6d^* 

3p 

0 

29426 

402 

32 

-155 

29705 



29840^, 29833^ 



1 

30659 

560 

10 

-164 

31065 

31249 

-184 

31365^, 31342_^ 



2 

31892 

1257 

22 

-142 

33029 

32941 

88 

33180^, 33147^ 

7s8p 

1 po 

1 

30691 

51 

55 

-94 

30703 

32858 

-2155 

30695^* 

&d^ 

Ij) 

2 

30490 

493 

13 

-144 

30852 



30982^, 30930^ 

6d7p 

ipo 

2 

30747 

312 

43 

-150 

30952 

30918 

34 

30894^, 31178^ 

7s8p 

3pO 

0 

31180 

-152 

79 

-53 

31054 

31086 

-32 

31008^, 31126^ 



1 

31511 

-5 

55 

-88 

31473 

31563 

-90 

31446^*, 31636^ 



2 

31913 

-92 

75 

-64 

31832 

31874 

-42 

31778^, 31934^ 

7s7d 


1 

32079 

-144 

74 

-51 

31958 

32001 

-43 

31895^, 32423^ 



2 

32749 

-810 

69 

-65 

31943 

31993 

-50 

31902^ 



3 

32273 

-143 

77 

-51 

32156 

32197 

-41 

32068^, 32625^ 

7s7d* 

iR) 

2 

32171 

26 

53 

-84 

32166 

32215 

-49 

32205^,32564_^ 

6d7p 

3p° 

1 

32180 

-1 

52 

55 

32136 

32230 

-94 

32090^, 32614^ 



2 

32287 

279 

17 

-157 

32426 

32507 

-81 

32436^, 32846^ 



3 

33021 

254 

28 

-147 

33156 

33197 

-41 

33169^, 3353U 

6d7p 

CO 

0 

0 

33543 

247 

21 

-151 

33660 

33782 

-122 

33809^, 34055^ 



1 

33586 

239 

24 

-147 

33702 

33824 

-122 

33837^, 34102_^ 



2 

34193 

283 

43 

-147 

34372 

34383 

-11 

344211, 346771 

Qd^ 

If? 

4 

33142 

724 

-16 

-166 

33684 



332611 

Qd^ 

Is 

0 

33561 

312 

54 

-116 

33811 



339611 

6d7p 

Ipo 

3 

34230 

155 

55 

-121 

34319 



34332I 

7s9s 

®s 

1 

34601 

-224 

81 

-41 

34417 

34476 

-59 



NIST data, Refj42l|. 

^ CP+CI, Ref. I29I]. I Pl found in combination with the calculation for the ion. Ref. [4,'ill . 

^ XIHFSCC, Ref. I31II . The value for IPl includes a QED radiative and frequency-dependent Breit correction of —46cm“^ from Ref. 

Semi-empirical CP-I-CI, Ref. [^ . 

« IHFSCC, Ref. [Ji. 
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Ref. [ 2 ^. In Ref. [3l|, the value for IPl was modified 
by adding the frequency-dependent Breit and radiative 
QED corrections from Ref. [s^; the contribution from 
the radiative corrections, through use of the radiative po¬ 
tential [ 3 ^, is —42cm“^, in agreement with the current 
work. In Ref. [ 2 ^, while Breit and QED radiative cor¬ 
rections were not included explicitly, these effects were 
taken into account to some degree through the empirical 
fitting factors. 

The current calculations for radium are the most com¬ 
plete and most accurate to date. With the high accuracy 
that we have achieved using this method, we can be con¬ 
fident of resolving anomalies with experiment in favour 
of the theoretical predictions. 

The most striking disagreement is with the energy as¬ 
signed to the state 7s8p ^Pi- In the experimental data 
[4^ . the large excitation energy 32858cm“^ is given to 
this state. However, we do not see this level in our cal¬ 
culations; rather, we see a low-lying state, with energy 
30703 cm“^, that is absent in the data. In the relativis¬ 
tic regime, the states are defined by their total angular 
momentum J and their parity. We see only six levels 
with J = 1 and odd parity in the energy range we have 
considered, and even extending this range, the next level 
that we see appears as high as 36067 cm“^. 

As for other supporting theoretical calculations for the 
anomaly, results of the semi-empirical CP-I-CI calcu¬ 
lation [23 included two levels with J = 1 and odd parity 
with energies 30695 cm“^ and 31446 cm“^. These were 
assigned in that work to the two levels 7s8p ^P° and 
7s8p ^Pi in the experimental data, with very large devi¬ 
ations —868cm“^ and —1412cm“^ being seen there. We 
suggest that these calculated excitation energies should 
instead be assigned to the low-lying levels 7s8p ^Pi and 
7s8p ^Pi- These numbers are in very close agreement 
with the results of the current work, 30703 cm“^ and 
31473 cm“^. No data for the corresponding singlet state 
^P° was presented in the IHPSCC work [iq . 

Therefore, we suggest that there is no high-lying singlet 
state ^P° with energy 32858 cm“^, though we expect that 
there should be a lower-lying one with energy around 
30700 cm“^. 

We also suggest that configuration assignments for 
some terms be altered. We do not have any terms in 
our data with leading configuration 7p^ as appears in 
the experimental data. Por the triplet terms 7p^ ^P, our 
results indicate instead a strongly dominating configura¬ 
tion 6(P ^P. There is also a case where we see that the 
dominant configuration is 7s7d rather than 7p^, for the 
term referred in the experimental data as 7p^ ^D 2 . 

In a system as heavy as radium, and particularly for 
the higher levels, the validity of the LS system for desig¬ 
nating the terms loses much of its meaning. Nevertheless, 
we still believe that making the previous observations is 
important, especially when different assignment of the 
conhgurations may lead to confusion between the desig¬ 
nations of levels. 

We note further, that with the very strong mixing of 


conhgurations in some terms, these designations become 
less clear, as we have seen for the odd-parity J = 1 terms 
7s8p ^P°, 7s8p and 6d7p ^D°, where the mixing 
between 7s8p and 6d7p is very strong. 


3. Element 120 

Our results for the ionization potentials and excita¬ 
tion energies for E120 are presented in Table I VIII We 
also present the results of our calculations for g-factors, 
and the corresponding non-relativistic values, to help in 
identihcation of the levels. Por the higher levels, there 
is strong mixing between conhgurations, and we have in¬ 
cluded the dominant conhgurations explicitly in the ta¬ 
ble. 

It is seen that the ladder diagrams for E120 are sig- 
nihcant, around 400 cm“^ for many levels, although the 
maximum correction for the levels for E120 is less than 
we saw for the lighter atoms. The d orbitals are well- 
screened by the relativistically contracted s and pi /2 or¬ 
bitals, and there are no low-lying terms with conhgu- 
ration. 

The Breit corrections to the excitation energies are 
about 150 cm“^ for many levels, while the QED radiative 
corrections are smaller, reaching 158 cm“^ in the largest 
case, for the level 8p^ 

In the hnal columns we list the results of other ab ini¬ 
tio calculations: CP-bCI [ 2 ^, a version of PSCC [30l| . 
and XIHPSCC [3l|. We are also aware of another calcu¬ 
lation [ 3 ^ of the excitation spectra of E120, though the 
results are so different from all other data (by as much 
as 10000 cm“^ for some levels) that we have decided not 
to present them in the table. 

Overall, there is good agreement between the different 
calculations presented in the table. It should be noted 
that Breit and QED corrections were not included in the 
calculations [ 2 ^, while QED was not included in Refs. 
Sill. Also, there is some sensitivity to the choice 
of nuclear density. We used the same Fermi distribution 
that was used in the previous calculation [ 2 ^ , and in that 
work the volume isotopic shifts were studied; we refer the 
reader to that work for more details. 

The QED radiative correction to the ionization poten¬ 
tial IPl for E120 was calculated rece ntly . Shabaev et al. 
used their model operator approach to calculate the 
self-energy correction; they obtained — 202cm“^. The 
total QED correction (including vacuum polarization) 
was calculated in Ref. using the radiative potential 
method [s^, and they obtained a result of —77 cm“^, in 
excellent agreement with our own —75 cm“^. Their result 
is broken down into the self-energy and vacuum polariza¬ 
tion contributions, —183cm“^ and -1-106cm“^, with the 
former in good agreement with the result of Shabaev et 
al 

The only other calculation of QED corrections to exci¬ 
tation energies for E120 was performed in the work [ 3 ^. 
The corrections were found using the radiative potential 




TABLE VII. Ionization potentials (removal of one electron, IPl, and both electrons, IP1+IP2) and excitation energies for 
E120 up to 40000 cm~^. CP“*+SD+CI excludes ladder diagrams, Breit, and QED; these corrections are given in subsequent 
columns. Results of other calculations are also presented. Units: cm“'^. Calculated and non-relativistic g-factors are given in 
the final columns. 


Conf. 

Term 

J 

Cpnl 

-bSD-bCI 

Ladder 

Breit 

QED 

Total 

CP-bCP 

FSCC^ 

XIHFSCC^ 

g 

Snr 

TPI 


0 

47702 

-510 

565 

-75 

47682 

47356 

47046 

47089 



IPl-bIP2 Ss^ 

is 

0 

137807 

-1030 

457 

-203 

137031 

137501 

136332 

136920 



8s8p 

3 

0 

16039 

-136 

123 

-90 

15936 

15777 

15328 

15648 

0.000 

0.000 



1 

18035 

-148 

125 

-92 

17920 

17710 

17382 

17587 

1.409 

1.500 



2 

25924 

-338 

153 

-88 

25651 

25419 

25308 

25192 

1.499 

1.500 

8s7d 


1 

22904 

69 

28 

-136 

22865 

22985 

22337 

22903 

0.500 

0.500 



2 

23299 

76 

42 

-134 

23283 

23163 

22494 

23034 

1.162 

1.167 



3 

23763 

96 

79 

-125 

23813 

23799 

23377 

23782 

1.333 

1.333 

8s8p 

Ipo 

1 

27827 

-331 

142 

-79 

27559 

27667 

28304 

27513 

1.087 

1.000 

8s7d 

ip) 

2 

27666 

41 

95 

-133 

27669 

27438 

27652 

27247 

1.006 

1.000 

8s9s 

®S 

1 

31464 

-427 

145 

-63 

31119 



30862 

1.998 

2.000 

8s9s 

Is 

0 

32633 

-431 

148 

-64 

32286 




0.000 

0.000 

8s9p 

3pO 

0 

35974 

-428 

159 

-75 

35630 



35463 

0.000 

0.000 



1 

36109 

-427 

159 

-77 

35764 



35595 

1.369 

1.500 



2 

37908 

-462 

163 

-76 

37533 



37369 

1.500 

1.500 

8s8d 


1 

37880 

-427 

152 

-82 

37523 




0.501 

0.500 



2 

37932 

-426 

151 

-81 

37576 




1.145 

1.167 



3 

38094 

-424 

157 

-79 

37748 




1.333 

1.333 

8s9p 

1 po 

1 

38131 

-457 

161 

-77 

37758 




1.129 

1.000 

8s8d 

ip) 

2 

38341 

-417 

153 

-85 

37992 




1.022 

1.000 

8sl0s 

®s 

1 

39441 

-507 

159 

-69 

39024 




1.999 

2.000 

8 pV8s10s 

Is 

0 

39448 

-374 

180 

-112 

39142 




0.000 

0.000 

8 pV8s10s 

Is 

0 

40137 

-341 

206 

-158 

39844 




0.000 

0.000 

8s7f/8s6f 

3p° 

2 

40139 

-296 

167 

-104 

39906 




0.680 

0.667 



3 

40410 

-352 

168 

-77 

40149 




1.061 

1.083 



4 

40407 

-352 

168 

-76 

40147 




1.250 

1.250 

8s7f/8s6f 

IpO 

3 

40255 

-333 

170 

-83 

40009 




1.024 

1.000 


Ref. IPl found in combination with the calculation for the ion, Ref. [4.111 . 

^ Ref. M. 

Ref. The value for IPl includes a QED radiative and frequency-dependent Breit correction of —102cm“^ from Ref. 1371 


method implemented in GRASP [13, IH, H^, with 
results in very good agreement with the results of this 
work. 


4- Accuracy 

The quality of our calculations can be gauged by com¬ 
parison of our results with experiment. We found excel¬ 
lent agreement between theory and experiment for bar¬ 
ium and radium, 100 cm“^ for many cases, and we can 
expect a similar level in the error for our predictions for 
E120. 

We note that the largest deviations we have seen be¬ 
tween theory and experiment mostly involve configura¬ 
tions containing d-orbitals. The smaller deviations we 
see for Ra compared to Ba could be due to the d-orbitals 
being screened (also there are fewer configurations involv¬ 
ing d-orbitals in the lower levels); this screening effect is 
even more pronounced for E120, and therefore the asso¬ 
ciated errors from d-orbitals may be reduced. 


To further study the limits of our method, particularly 
in relation to the quality of the correlation potential, we 
have also performed calculations of the spectra for the 
three atoms using a different set of screening factors for 
the Coulomb interaction in the exchange part of the cor¬ 
relation potential. The screening factors were found from 
the direct diagrams with both electron-electron screening 
and the hole-particle interaction included, see the sec¬ 
ond row of Table H] Calculations performed with these 
modified screening factors gave results that differed from 
our final results “Total” in Tables |V1 IVIl IVIII over the 
full spectral range we considered by about 10 — 200 cm“^ 
for barium, about 10 — 100 cm“^ for radium, and about 
10 — 200 cm“^ for element 120. For E120, the largest 
deviations occurred for the higher levels corresponding 
to configurations comprised of valence orbitals that are 
not the lowest for that wave. For many of the levels for 
Ba and Ra, the results with the adjusted screening fac¬ 
tors were better than or as good as the original results 
(and often on the other side of the experimental values). 
Therefore, we estimate an error that is about half the size 
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of these shifts, consistent with the error estimates we ob¬ 
tained by looking at the deviation of the original results 
from experiment. We expect similar behaviour for E120. 

IV. CONCLUSION 

In summary, we have performed ab initio calculations 
of the spectra of Ba, Ra, and E120 using the recently- 
developed CP-I-SD-I-CI method. We have found unprece¬ 
dented agreement with experiment for Ba and Ra and 
have made accurate predictions for missing and unreli¬ 
able data and for the spectra of E120. For Ra, we are 
confident that the energy assigned to the level 


in the data is incorrect, and we predict a value 
Ri 2200 cm“^ smaller. The size of the error for our cal¬ 
culations for Ra is about 100 cm“^ or better, while it is 
slightly larger for Ba. We estimate a similar level of un¬ 
certainty for the spectra of E120. Finally, we note that 
the size of the errors in this ab initio method is compa¬ 
rable to the size of the QED radiative corrections. 


ACKNOWLEDGMENTS 

This work was supported in part by the Australian 
Research Council. 


[1] R. J. Holt et ai, Nucl. Phys. A 844, 53c (2010). 

[2] L. Willmann et ai, Hyp. Int. 211, 39 (2012). 

[3] V. Spevak, N. Auerbach, and V. V. Flambaum, Phys. 
Rev. C 56, 1357 (1997). 

[4] J. Dobaczewski and J. Engel, Phys. Rev. Lett. 94, 232502 
(2005). 

[5] W. C. Griffith et ai, Phys. Rev. Lett. 102, 101601 (2009). 

[6] V. V. Flambaum, Phys. Rev. A 60, R2611 (1999). 

[7] V. A. Dzuba, V. V. Flambaum, and J. S. M. Ginges, 
Phys. Rev. A 61, 062509 (2000). 

[8] J. S. M. Ginges and V. V. Flambaum, Phys. Rep. 397, 
63 (2004). 

[9] V. A. Dzuba, V. V. Flambaum, J. S. M. Ginges, and M. 
G. Kozlov, Phys. Rev. A 66, 012111 (2002). 

[10] J. Bierori et ai, Phys. Rev. A 80, 012513 (2009). 

[11] K. V. P. Latha and P. R. Amjith, Phys. Rev. A 87, 
022509 (2013); Phys. Rev. A 87, 039906 (2013). 

[12] B. M. Roberts, V. A. Dzuba, and V. V. Flambaum, Phys. 
Rev. A 89, 042509 (2014). 

[13] L. Radziute, G. Gaigalas, P. Jonsson, and J. Bieroh , 
Phys. Rev. A 90, 012528 (2014). 

[14] S. M. Ramachandran and K. V. P. Latha, Phys. Rev. A 
90, 042503 (2014). 

[15] J. Bieroh, G. Froese Fischer, S. Fritzsche, and K. 
Pachucki, J. Phys. B 37, L305 (2004). 

[16] J. Bieroh and P. Pyykko, Phys. Rev. A 71, 032502 (2005). 

[17] J. Bieroh, J. Phys. B 38, 2221 (2005). 

[18] A. Landau, E. Eliav, Y. Ishikawa, and U. Kaldor, J. 
Ghem. Phys. 113 9905 (2000). 

[19] V. A. Dzuba and J. S. M. Ginges, Phys. Rev. A 73, 
032503 (2006). 

[20] V. A. Dzuba and V. V. Flambaum, J. Phys. B 40, 227 
(2007). 

[21] J. Bieroh, P. Indelicato, and P. Jonsson, Eur. Phys. J. 
Spec. Top. 144, 75 (2007). 

[22] N. D. Scielzo et al, Phys. Rev. A 73, 010501 (R) (2006). 

[23] W. Trimble et ai, Phys. Rev. 80 054501 (2009). 

[24] B. Santra, U. Dammalapati, A. Groot, K. Jungmann, and 
L. Willmann, Phys. Rev. A 90, 040501 (2014). 

[25] V. A. Dzuba, Phys. Rev. A 90, 012517 (2014). 

[26] J. H. Hamilton, S. Hofmann, and Y. T. Oganessian, 
Annu. Rev. Nucl. Part. Sci. 63, 383 (2013). 

[27] K. Rutz et al, Phys. Rev. C 56, 238 (1997). 

[28] J. Even et ai, Science 345, 1491 (2014). 


[29] T. H. Dinh, V. A. Dzuba, V. V. Flambaum, and J. S. M. 
Ginges, Phys. Rev. A 78, 054501 (2008). 

[30] L. V. Skripnikov, N. S. Mosyagin, and A. V. Titov, Ghem. 
Phys. Lett. 555, 79 (2013). 

[31] A. Borschevsky, V. Pershina, E. Eliav, and U. Kaldor, 
Phys. Rev. A 87, 022502 (2013). 

[32] M. S. Safronova, M. G. Kozlov, W. R. Johnson, and D. 
Jiang, Phys. Rev. A 80, 012516 (2009). 

[33] V. A. Dzuba, V. V. Flambaum, and M. G. Kozlov, Phys. 
Rev. 54, 3948 (1996). 

[34] V. A. Dzuba, V. V. Flambaum, and O. P. Sushkov, Phys. 
Lett. A 140, 493 (1989). 

[35] V. A. Dzuba, Phys. Rev. A 78, 042502 (2008). 

[36] V. V. Flambaum and J. S. M. Ginges, Phys. Rev. A 72, 
052115 (2005). 

[37] C. Thierfelder and P. Schwerdtfeger, Phys. Rev. A 82, 
062503 (2010). 

[38] T. Hangele, M. Dolg, and P. Schwerdtfeger, J. Phys. 
Ghem. 138, 174113 (2013). 

[39] K. Dyall, I. Grant, C. Johnson, F. Parpia, and E. Plum¬ 
mer, Comput. Phys. Comm. 55, 425 (1989). 

[40] F. Parpia, C. F. Fischer, and I. Grant, Comput. Phys. 
Comm. 175, 745 (2006). 

[41] W. R. Johnson and J. Sapirstein, Phys. Rev. Lett. 57, 
1126 (1986). 

[42] A. Kramida, Yu. Ralchenko, J. Reader, and NIST ASD 
Team (2014). NIST Atomic Spectra Database (ver. 5.2), 
[Online]. Available: http://physics.nist.gov/asd Na¬ 
tional Institute of Standards and Technology, Gaithers¬ 
burg, MD. 

[43] T. H. Dinh, V. A. Dzuba, V. V. Flambaum, and J. S. M. 
Ginges, Phys. Rev. A 78, 022507 (2008). 

[44] V. A. Dzuba, Phys. Rev. A 88, 042516 (2013). 

[45] B. K. Mani and D. Angom, Phys. Rev. A 83, 012501 

( 2011 ). 

[46] H. P. Palenius, Phys. Lett. A 56, 451 (1976). 

[47] S. J. Rose, N. C. Pyper, and I. P. Grant, J. Phys. B 11 , 
755 (1978). 

[48] E. Eliav, U. Kaldor, and Y. Ishikawa, Phys. Rev. A 53, 
3050 (1996). 

[49] V. A. Dzuba and W. R. Johnson, Phys. Rev. A 57, 2459 
(1998). 

[50] V. M. Shabaev, 1. I. Tupitsyn, and V. A. Yerokhin, Phys. 
Rev. A 88, 012513 (2013). 



